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We combine tlie results of searclies by tlie CDF and DO collaborations for a fermiophobic Higgs 
boson {Hf) produced via the processes WHf, ZHf, and vector-boson fusion in pp collisions at 
-^s = 1.96 TeV at the Fermilab Tevatron collider. The analyses seek Higgs boson decays to W^W~ 
and 77. With up to 8.2 fb~^ of integrated luminosity analyzed at CDF and up to 8.2 fb~^ at DO, 
we obtain a 95% CL lower bound on the mass of the Higgs boson in the fermiophobic Higgs model 
of 119 GeV/c^ 
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I. INTRODUCTION 

The mechanism for the breaking of the SU (2) x U{1) electroweak gauge symmetry [l| of the standard model (SM) Q 
has thus far evaded experimental verification even though the SM accommodates current collider data and it has been 
tested with remarkable precision. Searches for the SM Higgs boson at LEP [3], the Tevatron [4], and the LHC [1, |i| 
have yet to confirm its existence. 

The mechanism of electroweak symmetry breaking may offer a richer phenomenology than expected in the minimal 
version of the SM. Several Higgs bosons may exist, or the Higgs boson(s) may have couplings that are different from 
those expected in the minimal model. Here we explore the possibility that the lightest Higgs scalar does not couple to 
fermions at tree level, but is in all other ways SM-like. In order to provide masses to the fermions, additional degrees 
of freedom must exist in the Higgs sector, and models with Higgs doublets and triplets are possible 04^- We consider 
here a model in which all particles and interactions beyond those of the SM have no phenomenological impact for 
observables at the Tevatron aside from suppressing the fermion couplings to the Higgs boson. We call this model the 
fermiophobic Higgs model (FHM). In the following, we denote the fermiophobic Higgs boson Hf, the SM Higgs boson 
HsM, and use the symbol _ff in a context which is valid for both. 

In the FHM, the production of Higgs bosons at hadron colliders via the process gg — > Hf is suppressed to a negligible 
rate. In the SM, this process proceeds through one-loop intermediaries via fermion couplings to the Higgs boson, and 
at two or more loops via fermion and boson couplings to the Higgs boson. Two-loop electroweak processes contribute 
to gg — Hf production in the fermiophobic Higgs model (FHM), but at a vanishingly small rate, as their main 
contribution in the SM are through interference with one-loop processes [lOj . The associated production mechanisms 
pp — >■ WHf + X and pp ZHf + X, as well as the vector-boson-fusion (VBF) processes q'q Hfq"q"' remain 
nearly unchanged relative to the corresponding processes in the SM. Table |T] lists the production rates we use for the 
WHf, ZHf, and VBF modes. We neglect the cross section for gg — > Hf in this model. We use the WHa,M and ZHsm 
cross sections from Ref. which are based on the next-to-leading-order (NLO) calculation of v2hv [l^l and include 
next-to- next-to-leading-order (NNLO) QCD contributions [l^, as well as one- loop electroweak corrections We 
use the VBF cross section computed in the SM at NNLO in QCD in Ref. [l^. Electroweak corrections to the VBF 
production cross section are computed with the hawk program [l^, and are small and negative (2-3%) in the Higgs 
boson mass range considered here. We include these corrections in the VBF cross sections used in this combination 
of results from the Tevatron. We use the systematic uncertainties on the production cross section predictions from 
their respective publications, treating the WHf and ZHf production cross section uncertainties as fully correlated 
with each other, and uncorrelated with the VBF production cross section uncertainty. 

Another consequence of the FHM is to modify the decay branching ratios of the Higgs boson with respect to the SM 
predictions. Direct decays to fermions are forbidden; decays to fermions can proceed only via loops involving weak 
gauge bosons, and these are highly suppressed due to a spin flip required in the loop. The partial width of the Higgs 
boson decay to a pair of gluons, which in the SM is comparable to that of the decay to a pair of charm quarks, is 
negligibly small in the FHM. The remaining decays, to 77, W^W~ , Z7, and ZZ, account for nearly the entire decay 
width. For Higgs boson masses below twice the W boson mass, m/f^ ^ 2My^, the favored direct decay to W^W~ 
is reduced because one or both of the W bosons is far from resonance. Nonetheless, because the competing decay 
modes have small partial widths, the branching to W^W^ is still dominant, particularly for a heavier Higgs boson 
with ruHf > 120 GeV/c^. The branching fraction B{Hf — > 77) is greatly enhanced over B{Hsm 77) for all mn, 
and its clean signature and excellent mass resolution provide most of the search sensitivity for ruHj < 120 GeV/c^. 

Table |I] lists the decay branching fractions we assume in the searches described below, all of which were computed 
with HDECAY [l^. The dominant uncertainties on the SM Higgs branching fractions in the range 100 GeV/c^ < 
ruHsm < 200 GeV/c^ arise from uncertainties on the b quark mass rub, the charm quark mass rric, and the strong 
coupling constant as [isl . [l9| . Because these uncertainties affect only the decay modes in which a fermion (specifically, 
a quark) couples to the Higgs boson, they do not affect the Higgs boson's branching ratios in the FHM, and therefore 
they do not contribute to this analysis. 

The CDF and DO collaborations have searched for the SM Higgs boson in Hsm 77 decays [1^ These 
searches cannot be directly reinterpreted as constraints on the FHM by simply scaling the cross sections and branching 
ratios, because the kinematic distributions of the Higgs bosons, their decay products, and the particles produced in 
association with the Higgs boson differ between the FHM and the SM. These differences arise from the absence of 
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TABLE I: The production cross sections for WHf, ZHf, and VBF at NNLO Also listed are decay branching fractions (B) of 
the Higgs boson in the Fermiophobic Higgs model computed with hdecay [l^]. 
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gg — >■ Hf production in the FHM, and the increased fraction of the presence of an associated vector boson W or Z, or 
recoiling quark jets in the case of VBF. The transverse momentum (pt) spectrum of the Higgs boson in the FHM is 
thus much harder than it is in the SM, affecting the signal acceptance and the abihty to separate the signal from the 
backgrounds. The Hf — >■ 77 analyses have therefore been reoptimized for the FHM, taking advantage of the higher 
Pt of the Higgs boson 2ll, 23 . These analyses are updates of previous searches for the Higgs boson in the FHM using 
3.0 fb"i of data at CDF Q and 2.7 fb"^ of data at DO (H]. The CMS collaboration has also sought Higgs bosons 
decaying to photon pairs, and sets a lower limit on mHf in the FHM of 112 GeV/c^ at the 95% CL [25||. ATLAS has 
sought the SM Higgs boson in the decay Hgu 77 in 1.1 fb"i of data d]. The four LEP coUaborations, ALEPH, 
DELPHI, L3, and OPAL, have searched for e+e" — > ZHf — > Z77, and from their combined results report a lower 
limit on ruHf of 108.2 GeV/c^ in the FHM at the 95% CL ^26]. 

In addition to the searches for Hf 77, we also combine CDF's and DCs searches for H W^W^ , taking 
advantage of its enhanced branching fraction in the FHM. As for H 77, the searches for Hsm ~^ W^W~ cannot 
be interpted directly in the FHM due to the different mixture of production modes. 

The CDF collaboration's i/sM W^W^ analyses [l^l keep separate account of the predictions from gg — > i?sM, 
W^^SM, ZHsM, and VBF in each of the contributing channels. The distributions of the neural-network discriminants 
of these analyses are reinterpreted in the FHM by setting the gg — > Hf component to zero and by scaling the remaining 
signal components by the ratio of branching ratio predictions B{Hf — >• W'^W~)/B{Hsm W'^W^). CDF's searches 
for opposite-charge dilepton events in the final state Hsm — >■ W~^W~ i~^i'i£'~D£/ (where £,£' = e, fi), are separated 
into categories based on the number of reconstructed jets accompanying the leptons and missing transverse energy in 
the event. The gg ^ H process, which contributes negligibly in the FHM, most frequently produces a Higgs boson 
without additional jets, while the WH, ZH, and VBF production modes more commonly produce Higgs bosons 
with additional jets. This division of events into jet categories optimizes the search in the FHM and does not require 
developing a separate set of analysis channels, unlike the case of Hf — >■ 77. The opposite-charge low-m^^ channel and 
the channels seeking H — >■ W'^W~, in which one W decays into a t lepton which then decays to hadrons-|-^'T, while 
the other W decays to or fj,v, are included in the combination. The trilepton and like-charge dilepton analyses. 
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TABLE II: Luminosity, explored mass range and references for the different processes and final states = e or n) for the CDF 
analyses. The label "2x" refers to a separation based on lepton categories. 



Channel 


Luminosity mHf range 


Reference 




(fb-i) 


(GeV/c^) 




CDF Hf "f"f 


7.0 


100-150 


[21] 


CDF H W+W 2x(0 jets, I jet) + (2 or more jets) + (low-m«) + (e-rhad) + (/i-Thad) 


8.2 


110-200 


[27] 


CDF WH WW^W (same-sign leptons) + (tri-leptons) 


8.2 


110-200 


[27] 


CDF ZH -5> ZW'*'W~ (tri-leptons with 1 jet)-|-(tri-leptons with 2 or more jets) 


8.2 


110-200 


[27] 


DO Hf — >• 77 


8.2 


100-150 


[22] 


DO VH e^e^ + X 


5.3 


115-200 


[28] 



also included in CDF's SM H — > W^W~ search [27], are included here as well, as they are targeted at the WH and 
ZH signal contributions. 

The DO SM searches for iJsM W'^W' included here are the like-charge dilepton searches [1^, targeting the WH 
and ZH production modes, in which the W or Z produced in association with the Higgs boson decays leptonically 
and forms a like-charge pair with a lepton from one of the W decays. They are reinterpteted here as searches in 
the FHM by scaling by the ratio of branching ratios in the SM and FHM. CDF uses neural-network discriminants 
and DO uses boosted decision trees as the final step in separating a possible signal from the backgrounds. The same 
multivariate discriminant functions used in the SM searches are used here without reoptimization for the FHM. A 
summary of the included channels and their luminosities is provided in Table HIl 

We use the same statistical methods employed in Ref. ^29], namely the modified frequentist (CLs) and Bayesian 
techniques, in order to evaluate the results. Pseudo-experiments drawn from the background predictions varied ac- 
cording to their systematic uncertainties are used to compute the limits we expect to obtain in the absence of signal. 
Correlated systematic uncertainties are treated in the same way as they are in Ref. [29| . The sources of correlated 
uncertainty between CDF and DO are the total inelastic pp cross section used in the luminosity measurement, the SM 
diboson background production cross sections {WW, WZ , and ZZ), and the tt and single top quark production cross 
sections. Instrumental effects such as trigger efficiencies, photon and lepton identification efficiencies and misidentifi- 
cation rates, and the jet energy scales used by CDF and DO remain uncorrelated. To minimize the degrading effects of 
systematics on the sensitivity of the search, the signal and background contributions are fit to data by maximizing a 
likelihood function over the systematic uncertainties for both the background-only and signal-|-background hypothe- 
ses [1^. We include the theoretical uncertainties on the predictions of the production cross sections when quoting 
limits normalized to them and when quoting our limit on tuhj in the FHM. 

The combined limits on Higgs boson production normalized to predictions of the FHM are listed in Table IHII for 
both the CLs and the Bayesian methods, for 100 GeV/c^ < ruHf < 200 GeV/c^, in 5 GeV/c^ steps. The expected and 
observed limits for both methods agree within 5% for nearly all rriHf values except for rriHf > 185 GeV/c^, where the 
agreement still remains close for the expected limit and is within 10% for the observed limits. We choose to present 
the limits calculated with the Bayesian approach which was selected a priori. The final limits are shown in Fig. [T] To 
quote a limit on in the FHM, we interpolate the limits linearly between the sampled values of ruHj and report 
the locations at which the observed and expected limit functions cross unity. We exclude niHf < 119 GeV/c^ at the 
95% CL, and the expected exclusion range is also niHf < 119 GeV/c^. 

In summary, we present a combination of CDF and DO searches for the Higgs boson in a model in which the 
tree-level couplings of the Higgs boson to fermions vanish, while retaining standard model couplings to bosons. We 
exclude, at the 95% CL, a Higgs boson of mass niHf < 119 GeV/c^ in this model. This is the most restrictive limit 
to date on the fermiophobic Higgs model. 
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TABLE III: Observed and expected limits from the combined Tevatron search for Higgs boson production in the fermiophobic 
Higgs model (FHM) 
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